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Abstract 



We investigate finite volume effects in the propagators of Landau gauge Yang- Mills 
theory using Dyson- Schwinger equations on a 4-dimensional torus. In particular, 
we demonstrate explicitly how the solutions for the gluon and the ghost propaga- 
tor tend towards their respective infinite volume forms in the corresponding limit. 
This solves an important open problem of previous studies where the infinite vol- 
ume limit led to an apparent mismatch, especially of the infrared behaviour, be- 
tween torus extrapolations and the existing infinite volume solutions obtained in 
4-dimensional Euclidean space-time. However, the correct infinite volume limit is 
approached rather slowly. The typical scales necessary to see the onset of the lead- 
ing infrared behaviour emerging already imply volumes of at least 10 to 15 fm in 
lengths. To reliably extract the infrared exponents of the infinite volume solutions 
requires even much larger ones. While the volumes in the Monte-Carlo simulations 
available at present are far too small to facilitate that, we obtain a good qualitative 
agreement of our torus solutions with recent lattice data in comparable volumes. 
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1 Introduction 



The infrared behaviour of the QCD Green's functions is known to contain essential information 
about the realisation of confinement in the covariant continuum formulation of QCD in terms 
of local field systems [1]. To extract this information from studies of correlation functions in 
a finite volume it is most important to address the volume dependence of the long-range, i.e. 
infrared, behaviour of these correlations. It is only when this dependence is under control that 
firm conclusions can be drawn from infinite volume extrapolations. 

In this paper we revisit the elementary 2-point correlation functions, i.e. the propagators, of 
pure Yang-Mills theory without quarks in a finite volume. Using covariant gauges, and here 
in particular the Landau gauge, these are the gluon and the ghost propagators. These are 
perhaps the most important examples of how Green's functions relate to confinement in this 
formulation. Their infrared behaviour is sensitive to confinement according to the scenarios of 
Kugo and Ojima [2,3], and of Gribov and Zwanziger [4,5]. In Landau gauge these scenarios 
predict an infrared vanishing gluon propagator and an infrared enhanced ghost propagator (the 
latter was pointed out in [6]; for a review, see [1]). 

Qualitatively, finite-volume effects are expected to distort the infrared properties of 2-point 
correlation functions for momenta approaching 2n/L, where L is the finite length of the system 
in the corresponding direction. For pure Yang-Mills theory the momentum scale for confinement 
is set by Aq C d, which is of the order of 200 - 400 MeV, somewhat dependent on the scheme0 
This implies that for lengths around L c = 27t/Aqcd, of the order 3-6 fm, the finite volume 
will already have a considerable effect on the infrared behaviour of these confinement-sensitive 
correlations. Eventually, of course, when further decreasing the volume down to the sub-Fermi 
regime, confinement will be lost altogether. 

In present Monte-Carlo simulations, using lattice Landau gauge, the gluon and ghost propaga- 
tors of pure SU(2) and £77(3) gauge theory in four dimensions are obtained in volumes of the 
order of L c in length, in some cases even somewhat larger than L c [8,9,10,11]. However, despite 
clear indications in favour of the qualitative infrared behaviour of both the gluon and the ghost 
propagator, as predicted by continuum studies, some quantitative discrepancies still remain. In 
particular, apart from some indications from the volume scaling of the zero momentum gluon 
propagator, a systematic and statistically significant verification of its infrared suppression from 
lattice data in 4 dimensions is still lacking] 2 1 

As far as the continuum studies are concerned, on the other hand, we are nowadays in the quite 
comfortable situation that a variety of different non-perturbative approaches all lead to the 
same infrared behaviour for the propagators in the infinite volume limit. These include studies 
of their Dyson-Schwinger Equations (DSEs) [14], of the Fokker-Planck type diffusion equa- 

1 E.g., Ajjs 240 MeV for Nf = from the lattice determination by the ALPHA Collaboration [7]. 

2 This is different in Coulomb gauge and interpolating gauges, where a stronger suppression is expected 
from the Gribov- Zwanziger scenario, and where this infrared suppression of the gluon propagator is 
observed on the lattice already in rather small volumes [12,13]. 
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tions of Stochastic Quantisation [15], and of the Functional Renormalisation Group Equations 
(FRGEs) [16]. 

In order to understand the origin of the remaining discrepancies observed between the functional 
methods in the continuum and the lattice Landau gauge simulations, it is an obvious and 
necessary step forward to adapt the continuum methods to finite volumes. The techniques 
to solve Dyson-Schwinger equations on a finite four dimensional torus with periodic boundary 
conditions have recently been developed and applied to the various propagator DSEs of Landau 
gauge QCD [17,18,19]. The infinite volume limit, however, remained unclear in these studies. 
Here we specifically address this question and explicitly demonstrate how the infrared behaviour 
known from the previous studies of functional methods is approached in the infinite volume 
limit also by the DSE solutions on a finite torus. 

In fact, this is to be expected. The perturbative and the confining regime are separated by a 
cross-over, the scale of which is Aqcd, the dynamical scale. As pointed out in previous infrared 
studies [14,16,20], for the investigation of the infrared regime of the correlation functions we 
have to consider Aqcd an ultraviolet scale which is large compared to all external momenta p 
involved, i.e., p <C Aq CD . As mentioned above, for the emergence of continuum infrared physics 
we need in addition that 2n/L <C p. To reliably extract infrared critical exponents from finite 
volume studies we therefore need a clear separation of scales with a sufficient number of different 
momentum values all in the range 

2tt 

— <p<A QCD . (1.1) 



Indeed, the same scale separation was used in the infrared analysis of the functional RGEs 
[16,20], where an explicit infrared momentum cut-off k is introduced and plays a role analogous 
to that by the finite volume. In the functional RG infrared-studies it is precisely this regime, 
k <^p <^ Aqcd, that is required to lead to the same infrared behaviour as in the other functional 
continuum investigations without such infrared cut-off [14,15]. In the available studies of the 
DSEs on a symmetric hypertorus, on the other hand, this same limit has not been explicitly 
obtained as yet. In the present work we fill this gap by presenting a refined infrared analysis. 

Because of the necessary separation of scales (1.1), the actual lattice sizes L necessary to achieve 
this separation for a sufficiently wide range of momentum values need to be much larger than 
L c . From our numerical analysis we estimate that this might well require lengths L of the order 
of 40 fm, so that reliable lattice determinations of the infrared exponents of the gluon and ghost 
propagators on reasonably fine lattices are likely to remain unfeasible for some time to come. 
Corresponding 577(2) simulations on rather large lattices in three dimensions [21] confirm this 
by showing a slow but steady convergence towards the three dimensional infinite volume limit, 
as predicted by DSE studies [15,22]. 

Meanwhile, we explicitly demonstrate this convergence here for the DSE solutions in a finite 
4-dimensional volume by investigating the behaviour of the propagators on length scales in the 
range given by (1.1). Sec. 2 is devoted to analytic results. After a brief summary of the infrared 
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behaviour of the propagators from infinite volume studies in Sec. 2.1, we discuss the torus 
DSEs and their proper renormalisation in Sec. 2.2. The question of zero modes is addressed 
in Sec. 2.3. The infrared behaviour of the torus-DSE solutions and the volume dependence 
of the low-momentum modes are studied in Sec. 2.4. Our numerical results are presented in 
Sec. 3. This contains a brief description of the general procedure and the implementation of 
the renormalisation procedure in Sees. 3.1 and 3.2, and a detailed discussion of the numerical 
solutions for the gluon and ghost propagators in various volumes in Sec. 3.3. In Sec. 3.4 we 
present fits and extrapolations to model the approach towards the infinite volume limit, and in 
Sec. 3.5 we compare our solutions to the results of lattice simulations with similar volumes. 

Our conclusions are provided in Sec. 4. We summarise that the presently available lattice data is 
consistent with the continuum results, if the consequences of finite-volume effects are properly 
taken into account. Moreover, the infinite-volume extrapolations of our torus-DSE solutions 
indicate that volumes of lengths of the order of 40 fm might be necessary to explicitly verify 
the infrared behaviour of the gluon and ghost propagators in the various functional continuum 
approaches, and to obtain reliable estimates from lattice simulations of the critical exponents 
indicative of this conformal infrared behaviour of the pure gauge theory in the infinite volume, 
the 4-dimensional Euclidean space-time. 



2 Analytic Results 

2.1 Dyson-Schwinger equations at infinite volume 

The Dyson-Schwinger equations of QCD in 4-dimensional Euclidean space-time have been used 
very successfully to determine the Green's functions in particular of the Landau gauge [1,23]. 

For the propagators of the pure gauge theory in the Landau gauge the full momentum depen- 
dence is parametrised by two dimensionless dressing functions Z and G as follows, 



where D^ v is the gluon propagator and D is the propagator of the Faddeev-Popov ghosts. 

The infrared behaviour of these dressing functions Z and G for gluons and ghosts is determined 
by one unique infrared exponent < k < 1 [24,25], 




D ab (p 2 ) 



■ ab G(p) 





(2.2) 
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The exact value of this exponent k depends only on the infrared behaviour of a single invariant 
function A(k 2 ; p 2 , q 2 ) [14] which multiplies the tree-level structure in the full ghost-gluon vertex. 
Here k 2 denotes the gluon momentum andp 2 and q 2 are the (anti)ghost momenta. In ghost/anti- 
ghost symmetric gauges, such as the Landau gauge, A is symmetric under the exchange of the 
ghost momenta, A(k 2 ; q 2 ,p 2 ) = A(k 2 ;p 2 , q 2 ). If this function is furthermore assumed to be finite 
and regular in the origin ctt k = p 2 = q 2 = 0, then the value of the infrared exponent k is [14] 



93 - VmK . . 

k = k c = « 0.595 . 2.3 

98 v ' 

This regularity assumption is all that is needed to obtain this otherwise exact result. The same 
value of k was obtained independently at the same time using a bare vertex (which is trivially 
regular) from the time-independent (equilibrium version of the) diffusion equation of Stochastic 
Quantisation [15]. In the corresponding truncation of the functional RGEs a range of values 
roughly between 0.54 and k c is possible with (2.3) representing the special value obtained for 
an optimised flow [16,20]. 

The same regularity assumption implicitly underlies all these studies and produces the unique 
result (2.3). Other values for k in the range 1/2 < k < 1 are possible, if the tree-level structure 
A of the ghost-gluon vertex has infrared divergences associated with any of its legs [14]. 

So how justified is this regularity assumption? It is relatively easy to see that A(p 2 ;0,p 2 ) = 
A(p 2 ;p 2 ,0) = 1 in Landau gauge. An independent and different argument, from the non- 
renormalisation of the ghost-gluon vertex in Landau gauge [26], perturbatively at all orders, 
implies A(p 2 ;p 2 ,p 2 ) = 1 in a symmetric momentum subtraction scheme. If this argument re- 
mains true beyond perturbation theory it fixes the infrared limit of A along a second direction. 
Lattice simulations of the ghost-gluon vertex have so far tested the line A(0;p 2 ,p 2 ) in £77(2) 
and SU (3) for momenta p down to values of the order of Aq CD , and obtained results consistent 
with A(0;p 2 ,p 2 ) = 1 with no systematic indications of any momentum dependence along this 
third direction either [27,28,29] 00 While none of this can prove the infrared regularity of this 
ghost-gluon vertex structure, there is certainly no evidence to the contrary at present. 

Probably more importantly, however, there is a self-consistent scheme to solve the full hierarchy 
of DSEs in the pure gauge theory asymptotically in the infrared [31]. Remember that the 
complete system of DSEs forms an infinite set of coupled non-linear integral equations between 
all n-point Green's functions. The self-consistent infrared-asymptotic solution to the whole 
tower of n-point Green's functions obtained in [31] reflects the infrared fixed-point behaviour of 
the pure gauge theory. The ghost-gluon vertex is regular in the infrared in this solution, which 
incorporates (2.2) and leads to a simple counting scheme for generalised power laws of this kind 
for all n-point Green's functions. Moreover, it has recently been proven that this solution is 
unique among this class of parametrisations of a conformal infrared behaviour [32]. 

Because of the non-renormalisation of the ghost-gluon vertex in Landau gauge, the ghost and 

3 For more results in SU(2), including some different kinematic regimes, see [30]. 
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Fig. 1. The propagator DSEs. Curly lines are gluons and dotted lines are ghosts. Lines with a large 
filled dot represent full propagators. Vertices with a small dot are bare and with a large dot are full 
and constructed in this truncation. 

gluon dressing functions provide a non-perturbative definition of the running coupling in this 
gauge via the RG invariant product [25] 



(2.4) 



where /i 2 is the renormalisation point. It approaches a finite positive value for p 2 — > in four 
dimensions for all < k < 1. The value of this infrared fixed-point depends on k, i.e., on the 
infrared properties of the ghost-gluon vertex (it vanishes at both end points of this interval, for 
k — > and for k — > 1). Moreover, the maximal value of the coupling at the infrared fixed-point 
is that obtained for the infrared-regular ghost-gluon vertex, with k = k c from (2.3), which 
yields [14], 



8tt T 2 (k c - 1)T(4-2k c ) 4tt . , 

for N c colours. This implies a c ~ 4.46 for SU{2) and a c ~ 2.97 for SU(3). 

Numerical solutions of the DSEs can only be obtained in specific truncation schemes. For the 
coupled system of ghost and gluon propagator DSEs, given diagrammatically in Fig. 1, one such 
scheme has been denned in [17,34]. It uses a bare ghost-gluon vertex, a choice well justified by 
both lattice studies [27,28,29,30] as well as calculations using the vertex DSE [33]. Furthermore 
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an Ansatz for the three-gluon vertex has been used that ensures the correct infrared and 
ultraviolet limits of the solutions for the ghost and gluon propagators. Contributions involving 
the four-gluon interaction have been neglected, a choice which only affects the intermediate 
momentum regime [1,23]. The resulting numerical solutions of the ghost-gluon system have 
been reported in [34] and are included here again in Sec. 3. In the following we will denote 
these solutions as 'infinite volume solutions' as opposed to the 'torus solutions' discussed in 
the next sections. Note that for the latter we assume the continuum limit in a finite volume 
even though, strictly speaking, our numerical solutions in Sec. 3 are obtained on a discrete 
and finite set of momentum values. The continuum limit thereby corresponds to removing the 
ultraviolet cut-off, of course. In Sec. 3.3 we will verify that the residual cut-off dependences of the 
propagators are negligible. Unlike lattice simulations, where disentangling discretisation errors 
and finite volume effects can be a rather challenging task, it is therefore not necessary to study 
the continuum limit more carefully for the DSE solutions, provided a consistent renormalisation 
scheme is implemented. 



2.2 Dyson-Schwinger equations on a torus 



For the formulation of the Dyson-Schwinger equations of Fig. 1 in a 4-dimensional hypercubic 
volume L 4 of length L in all directions, with periodic boundary conditions, the momentum 
integrals in the infinite volume DSEs are replaced by sums, 

d 4 a 1 



where the four dimensional vector of integers n G Z labels the discrete momentum values q n = 
(2n/L) n. In the following we truncate the DSEs shown in Fig. 1 such that the diagrams with 
(bare) four-gluon vertices are neglected, analogous to the truncation scheme used in the infinite 
volume computations as mentioned above. Denoting the gluon and ghost dressing functions by 
Z L and Gl in the finite volume, their corresponding DSEs become, 



„ g 2 N c Z 1 l q n V(k)p n G L (q n ) G L (p n ) 
Z3 q 71 2^ 



Z L (k) 3 L 4 ^ k 2 q 2 n p 2 n 



g 2 N C Z X 1 H 3g (p n , q n , k) N(p n , q n , k) Z L (q n ) Z L (p n ) 

2^ h2 n 2„2 > \ lJ > 



3 L 4 „ k 2 q 2 n p 2 n 

1 7 7 1 X- kV (Pn)qn G L (q n )Z L ( Pn ) 

= Z 3 -g N C Z X — ^ , (2.8) 



G L (k) k 2 q 2 n p n 



with p n = k — q n . We furthermore use the translational invariance of the sums under q n — > 
q n + 2irm/L with m £ Z 4 . The abbreviation kV(p)q denotes a contraction with the transverse 
momentum tensor, i.e. kV(p)q = k^P^ v {p)q v with P^ u {p) = S^ u —p^p u /p 2 . The integration kernel 
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N(p, q, k) and the Ansatz for the dressing function H 3g (p, q, k) of the three-gluon vertex can be 
found in Refs. [17,18]. Note also that the above equations are not restricted to external momenta 
k n = (2ir/L)n, and k may or may not be one of these discrete values. For k ^ (2n/L)n, the 
Equations (2.7), (2.8) should be symmetrised in q n and p n . 

The renormalisation is done in the infinite volume limit, which suffices for finiteness, analogously 
to finite temperature field theory. This can be shown on the level of the functional DSEs by 
adding and subtracting the infinite volume DSEs and identifying the renormalisation constants 
Zi G {Zi, Zi, Z 3 , Z 3 } at infinite volume with those in a finite volume Z i L , 



Z t , L = Zi . (2.9) 

Such a scheme facilitates the approach towards infinite volume. We emphasise that the scheme 
does not imply the same RG conditions at finite and infinite volume. Applying the renormali- 
sation scheme (2.9) to the propagator DSEs (2.8) and (2.7), and furthermore using the result 
Zx — lva. Landau gauge [26], we obtain 



1 1 9 2 N C ( 1 ^q n V(k)p n G L (q n )G L ( Pn ) f d\ qV(k')p' G(q) G(p' 



E 



-/ 



Z L (k) Z(k' 2 ) 3 \L 4 V k 2 q 2 n p 2 n J (2tt) 4 k' 2 q 2 p' 2 

g 2 N c Z 1 ( 1 ^ H 3g (p n ,q n ,k) N(p n ,q n ,k) Z L (q n ) Z L (p n ) 



L 4 



E 



k 2 q 2 n p 2 n 



d 4 q H 3g (p>, q ,k>)N(p', q ,k')Z(q)Z(p>y 



(2tt)< 



k' 2 q 2 p 



,12 



(2.10) 



1 1 _ 2 N [}_ ST kV (Pn)q n G L (q n ) Z L (p n ) 

G L {k) G{k' 2 ) 9 C U 4 n k 2 q 2 n p 2 n 

r d 4 q k'V{p')q G{q) Z{p') 
~ 1 (2tt) 4 k' 2 q 2 p' 2 

The sum of the renormalisation terms in each of the two equations is, of course, momentum 
independent. The dependence on k' (with p' — k' — q) is thus an illusion. Any momentum value 
can be used equally. A particular and convenient choice below will be k! = k (and p' = p). 

The RG-scheme in (2. 10), (2. 11) facilitates our discussion of finite volume effects in the next 
section. The subtraction is 0(4)-symmetric and hence keeps the 0(4)-symmetry violations 
by the finite volume at a minimum. Moreover, this subtraction scheme is not restricted to 
the ghost and gluon propagators. It can be extended self-consistently to the DSEs of general 
Green's functions and general truncation schemes [20]. The implementation of this scheme in 
our numerical treatment of the propagator DSEs is discussed in Sec. 3.2. 



(2.11) 
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2.3 Zero modes 



An important subtlety concerns the ghost propagator and vertex functions. In the covariant 
gauges there are constant ghost zero modes which decouple completely from the theory. This 
is a consequence of the unfixed global gauge symmetry. The action does not depend on these 
constant ghost/anti-ghost modes, as it only couples to derivatives of the ghost and anti-ghost 
fields. As a consequence, the Faddeev-Popov operator can only be inverted for non-vanishing 
momenta p 7^ 0. The ghost propagator is thus only defined with the corresponding projection 
such that G L (q) has the property that 



G L (q n )=G L {q n ){l-6 n0 ). (2.12) 

Eq. (2.12) entails that the zero ghost-momentum terms are absent in the propagator DSEs. Note 
that in contrast the gluon zero modes do contribute to the loop sums as the gluonic vertex 
functions do depend on the constant gauge field. However, transversality is not well-defined 
at vanishing gluon momentum. Hence, we have to separate the zero-momentum contribution 
explicitly, writing, 

D L M = 5, u L 2 C L 5 n0 + P, u {p n ) (1 - 5 n0 ) . (2.13) 



n 



Here Cl is some dimensionless constant which represents the value of the gluon propagator at 
momentum zero in units of L 2 . In principle, this constant is to be determined self-consistently 
from the set of algebraic equations representing the propagator DSEs for the discrete momenta 
in the periodic L 4 box along with the remaining Z(p n ). Up to the factor 4/3 from the tensor 
structure, the value of the gluon propagator at momentum zero in the finite volume is expected 
to be of approximately the same magnitude as for the lowest non-zero modes p n . This is what 
is typically observed in lattice simulations also. For the parametrisation (2.13) it implies that 

Z L (2n/L) 

° L ~ ■ (2 ' 14) 



Eqs. (2.12) and (2.13) are important to obtain well-defined algebraic equations from the DSEs. 
These equations are closed among themselves only when the external momentum coincides with 
one of the discrete momenta k = k m = (2it/L)m for m G Z 4 . With Eqs. (2.12) and (2.13) we 
rewrite the gluon DSE (2.7) for Z(k m ) with m ^ n and p mn = k m — q n , and separately for Cl 
with m = 0, as 



Zlk~r Z ' + ¥F ^ %j^ n>2 G L MG L (p m ,,) + (2.15) 

AL\K m ) OIj n ^ 0y1n K ra ( lnPran 

^ g 2 N c 1 x , m 2 n 2 — (mn) 2 . . . . 
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CZ^^OM + - = ^T^E^ H- ... . ,2,6, 

We did not repeat the contributions from the 3-gluon loop explicitly here again. Note that in 
the gauge- invariantly regularised full theory gauge invariance implies that the DSE (2.16) for 
Cl must be ultraviolet finite, and of the order given by (2.14). Hence there is no counter-term 
in the equation for Cl- The term from the ghost loop alone, which is given explicitly here, is 
not finite, of course. Nonetheless, its apparent ultraviolet divergence, which has to be cancelled 
by those of the terms not given explicitly here, does not affect the infrared analysis, as will be 
shown below. 

With the separation (2.13), the ghost DSE (2.8) for m ^ becomes, 



1 _y 9 2 N C v k 2 m q 2 n -(k m q n ) 2 9 2 N C C L 
r 7% )- Z3 ~ 1^ W^a G L {q n )Z L { Pn ) ^p— G L {k m ) . (2.17) 



There is no DSE for m = in this case. In the momentum range (1.1) of interest here, we have 
(27r) 2 <C L 2 k 2 m . With dressing functions G(k m ), Z(k m ) of the order one, and C L of the order 
of magnitude given by (2.14), the gluonic zero mode contribution to the ghost DSE (2.17) is 
therefore suppressed by an explicit factor of (27r) 2 /(L 2 k 2 n ) as compared to the other terms in 
this momentum range. 

An important practical consequence of the above analysis is that we can safely drop the con- 
tributions of the gluonic zero mode in the DSEs. This is important for our numerical analysis, 
c./., Sec. 3. 



2.4 Infrared analysis 



We proceed with an infrared analysis following that of [14,16]. Due to the missing 0(4)- 
symmetry the dressing functions and Gi(fc) are functions of the individual components 

ki, and not of k 2 alone. When using an 0(4)-symmetric ultraviolet cut-off in (2. 7), (2. 8), however, 
we have approximate 0(4) symmetry for momenta k 2 L 2 > 1. Therefore, in order to simplify 
the analysis, we may neglect the residual 0(4)-symmetry violation considering an 0(4)-average 
with respect to the external momentum k in the propagator DSEs. 

For k — > 0, with k 7^ k rn = (2ir/L)m, the gluon equation (2.15) only depends on Gi(q n ) and 
G L {p n ) -> G L {q n ) (where p n = k - q n , c.f., Sec. 2.2, and thus p n -> q n ) with q 2 n > Atx/L 2 . 
Note that this momentum regime is below the lowest momentum obtainable in the numerical 
calculations in Sec. 3. 

The ghost dressing functions in the ghost loop of (2.15) therefore remain finite in this limit (for 
monotonically decreasing Giiqn) the upper bound will be Gl(2tt/L)). Hence, we deduce from 
the right hand side of the gluon equation that Z(k)/k 2 has a finite yet direction-dependent 
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limit for k — > 0. The angle dependence in this limit relates to the missing 0(4)-symmetry on 
the torus. After angular averaging we thus write 



Z L (k 2 ) = J%± Z L (k) oc (l + 0{k 2 /» 2 z )) , (2.18) 

for k — > 0. Eq. (2.18) displays the behaviour of a dressing function of a particle with a screening 
mass fiz(L). Inserting this behaviour into the ghost equation (2.17) we analogously deduce for 
the ghost dressing function that 



G L (k 2 ) = J ^ G L (k) oc 1 + 0(k 2 /^ 2 G ) , (2.19) 

for k — > 0. In contrast to the gluon propagator, the ghost propagator remains massless in the 
finite volume. Due to the ghost screening mass ^g{L) the ghost dressing function Gi{k) is 
infrared finite at finite L, however. The infinite- volume renormalisation scheme of Eqs. (2.9), 
(2.10) and (2.11) will guarantee that both screening masses /iz/g(L) — > in the infinite volume 
limit L — > oo. 

By construction, the averaged dressing functions only depend on the magnitude of momentum. 
We emphasise that this approximation does not change powers in k 2 . In particular, it will not 
affect the infrared exponents of momenta k ^> 2n/L. In the regime (1.1) we can expand the 
dressing functions Z L (k 2 ) and G L (k 2 ) about their leading physical infrared behaviour, 



2\ _ v { QCD 



Z L (k 2 ) = Z m 



A 2 \ K z L l, 
QCD \ 

A 2 ^ Ka 




ki+^j k 2 + (j? z v va Q cd'^v; ' 



G ^=H^fk) l 1+tfG te*^JJ' (2 - 20) 



with 



5Z(0, 0) = 5G(0, 0) = , (2.21) 

and possibly A;-dependent k>z l /g l but constant Zm, Gir. For k — > we arrive at the limits 
(2.18) and (2.19) with subleading powers depending on k 2 /fi 2 z ^ G . The finite volume corrections 
are suppressed with l/(k 2 L 2 ), and the confinement scale works as an UV cut-off leading to 
corrections of the order /c 2 /Aq CD . The corresponding terms are hidden in the corrections 5Z 
and 5G respectively. We could also have absorbed the screening 'masses' (jl z /g m to the definitions 
of the SZ/ <5G-corrections, but we prefer to keep them explicitly in this infrared Ansatz for later 
convenience. 



For Aq CD 3> k 2 3> /jl g / z we are left with the physical infrared behaviour, Z L (k 2 ) ~ (k 2 )' 
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Giik 2 ) ~ (k 2 ) k °l. In particular, for L — > oo the parametrisations (2.20) tend towards the 
infinite volume infrared forms, 

G L (k>) - G(fc 2 ) = G m (l + SG (xi^.o)) , (2-22) 

with the limits 

— > —2k , k Gl — > k , and Hz/g ^ for L — > oo . (2.23) 
This follows immediately from the renormalised DSEs (2.10), (2.11). 

We are also quantitatively interested in the approach towards the infinite volume solutions, and 
concentrate on momenta 



- 2 >0, (2.24) 

A QCD 

which is achieved formally by taking the limit Aq CD — > oo. The limit of asymptotically large 
volumes, on the other hand, is obtained from k 2 — > oo. In general, the infinite- volume DSEs will 
require renormalisation also for Aq CD — > oo. This is different from the standard perturbative 
UV-renormalisation in that contributions of momenta at about Aqcd are removed. Once the 
infinite volume DSEs are properly renormalised, however, the limit k 2 — > oo can be performed 
without further subtractions in the torus DSEs (2.7) and (2.8): when the right hand sides in 
the form of (2.10) and (2.11) are used, the results remain finite. 

We first insert the parametrisations (2.20) into the gluon Dyson-Schwinger equation (2.10). 
After angular averaging, we consider the leading contribution from the ghost loop in the mo- 
mentum regime (1.1) which becomes 



a 2 \ -*z L ,2 , .,2 n 2]\T r 12 



Zm \k 2 + fi 2 zj k 2 3 L 4 



y g n v{k)g n ( a^ cd / a 2 qcb 



+ 



Herein, labels the loop momentum closest to k, i.e., the minimum of p^. For momenta in 
the range (1.1) the corrections due to SZ, 5G can be neglected. Moreover, because the ghost 
loop dominates in this momentum regime, we can safely drop all the other contributions, not 
given explicitly here. The infinite volume renormalisation vanishes due to the value of k c > 0.5 
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in (2.3): the momentum integral in (2.10) is finite and agrees with 1/Z(k' 2 ). We immediately 
conclude that 



lim K Gh {q 2 ) >0.5, (2.26) 

to guarantee finiteness of (2.25). In the approximate expression (2.25), Aq CD /A; 2 and fc 2 ///|^ G 
are all assumed to be sufficiently large but still finite. When these ratios tend to infinity we 
recover one of the infinite volume DSE conditions for the critical exponent k c (2.3) and the 
infrared fixed point a c (2.5) in [14]. 

Note also that the momentum sums in (2.25) receive their main contribution for momenta at 
about k 2 . Upon rearranging (2.25) we obtain up to sub-leading terms, 

A / A 2 \ K Zt +2«Gr r) 

4vr _ / Aqcd \ . , . 2tt 



Iz L (KG L ,k,—) , (2.27) 



a c N c \k 2 + pL 2 z J ^' ' L 

where 

2tt 167T 2 1 g TO P(% ra ( tf + ti Y ( tf + ti Y 

I Z L {K,K,—)-— T 2^ p - 2 „2 7 ..2 I „2 7 ..2 I ' 

and a c = (# 2 /47r) Zi R G^ R . Recall that Eq. (2.27) holds for < A; 2 < Aq CD . For its left hand 
side to have a non-vanishing and finite value in the limit (2.24) with Aq CD — > oo, we need to 
have at sufficiently large but still finite L, as in the infinite volume limit, 

k Zl + 2k Gl = . (2.29) 

Thus we define analogously, 

kg l = kl , k Zl = -2k l , (2.30) 

as in the L — > oo limit. Then, the dependence on the scale Aqcd is eliminated from the infrared 
analysis as necessary. We obtain, 

Air 9tt 

' I ZL ( KL ,k,^), (2.31) 



a c N c L 

and it is easy to verify that in the second limit, that of large volume via k 2 — > oo, this explicitly 
leads to the infrared condition (2.5) for a c from the infinite- volume gluon DSE in [14], as 

T , , 2 ^ t / \ 2 f d 4 q 1 qV(k)q fk 2 V L ( k 2 V L /o ooN 

^•^T^^^^/^^^b) y (2 - 32) 
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1 Y 2 {-k l )Y{2k l -1) 

2 T 2 (k l -1)T(4-2k l ) ' 

for k 2 — > oo. The condition in Eq. (2.31) therefore implies that k l — > n c in this limit, as expected. 

The infrared dominant contribution for Aq CD — > oo in the ghost DSE still needs ultravio- 
let renormalisation, as it does in the infinite volume case. With (2.20), (2.22) we obtain the 
corresponding contributions to the renormalised ghost DSE (2.11), with k' = k, 



A qcd ^ KL _ 1 (^qcd\ Kc 2 Nr 7 ( 1 v kV(p n )k 



Gm. \k 2 + &) G m \ k 2 ) y c in \L\^ nk k 2 q 2 n (p 2 n + ^z) 

r d\ kP(p)k /A|cdV c /A| 
7 (2vr) 4 fcVp 2 V 9 2 / 1 P 



A| CD V L f A Qcd V 2 ^ /■ A ( A IcbY C ( A lcA' 2Kc 



This equation is ultraviolet finite because the (weakly) momentum dependent exponent kl — > k c , 
whenever the momentum argument of the corresponding dressing function becomes large (of 
the order of A QCD ), as already derived from the ghost loop in the gluon DSE above. Eq. (2.33) 
can be rewritten in the form, 

4 9 / k 2 \ KL /A 2 \ k l- k c 

" hr— s- [I G (K c ,k,0)-I G (K c )} . (2.34) 



aN c ULK L \k 2 + fi 2 G J \ k 2 

Here, Ig(k, k, 2ir/L) is the corresponding loop sum (as explicitly given in (2.40) below), which 
is ultraviolet divergent in this case. This divergence is cancelled by the likewise ultraviolet 
divergent infinite volume integral expression 

whose finite contributions are independent of k and given by [14] 

t ( K \ 3 r 2 (-«)r(2«-i) 

IM = -2r(-2«)r(«-i)r(« + 3 ) " (2 ' 36) 

That leaves only the divergent part in the terms in brackets on the right of Eq. (2.34) to 
remove the ultraviolet divergence of the loop Ig(kl, k, 2n/L) (with k l —>k c at large momenta). 
Eq. (2.34) is the analogue to Eq. (2.31) from the gluon DSE. 

Alternatively, we can use k! = in the renormalisation terms of (2.11) to compute the infinite- 
volume ghost renormalisation constant Z 3 explicitly [14], 

9 3 f d A q Z(q 2 )G(q 2 ) 
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Using the infinite volume infrared forms (2.22) for asymptotically large Aqcd with 0(4)- 
symmetric ultraviolet cut-off Auv this becomes 



^ e (_£)"• (238) 

Of course, asymptotically large Aqcd in this case means Aqcd — ► Auv 5 and Aqcd acts as the 
ultraviolet momentum cut-off in the infrared analysis as explained above. Then, with (2.20) in 
the ghost DSE (2.8), we obtain from the leading infrared contribution in this limit 

lG L (KL,k,—) + -[ u , - 2 . (2.39) 



a c N c *w"™> L A\k 2 + ^ 2 G 

The loop sum herein is cut off at q 2 n < A QCD = {2n /L) 2 N\, 

2n _ IQn 2 " 2 |5 kV( Pn )k ( tf+j4 Y L hi + fi \ 



2k l 



(2.40) 



where for asymptotically large q 2 , approaching Aq CD , the sum over the loop-momentum q n 
tends to the corresponding integral because of the large number L A q 2 dq 2 /(16n 2 ) of modes in 
[<Zn> Qn + dq„]. The ultraviolet divergence of the loop with cut-off Aqcd is therefore subtracted 
correctly by the infinite volume counter-term in (2.39), which thus remains finite in the Aq CD — > 
oo limit as required. 

As for the ghost-loop contribution to the gluon DSE, c.f., Eq. (2.31) with (2.32), in the second 
limit of large L 2 k 2 we recover the infrared self-consistency condition from the infinite- volume 
ghost DSE of [14]. Using Eqs. (2.31) and (2.39), or (2.34), we furthermore obtain the leading 
corrections oc l/(k 2 L 2 ) as follows: 

First, expanding Iz L (K, L ,k,2n/L) in Eq. (2.28) for large k 2 , from the ghost loop of the gluon 
DSE (2.31), we obtain 



An 

Oic.Nc 



I z (n L ) |1 + (2« L - l)^ - 2n L {1 + Kl)(1 _ Kl) p ) + 0(j^) . (2.41, 



Analogously, from the ghost DSE, condition (2.39) becomes 



An 



Ig{k l ) 1 1 - [k l - 



k\ (2 + k l ) \ /4 
(1 + k l )(1-k l )J k 2 

(2ki 



(2.42) 



1)(2 + k l ) /4 

k 2 



4(1 - K L ) 
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The leading order herein again reproduces the corresponding infinite volume condition which 
entails kl — > n c in the limit k 2 — > oo, as 

i;T I z (k c ) = I g (k c ) . (2.43) 



a r .N, 



With J z (k), I g (k) from Eqs. (2.32), (2.36) it follows that k c = (93 - v / 120T)/98 = 0.59535... 
and a c = 2.9717... (for N c = 3), see [14]. 

Secondly, for the leading corrections to this limit, we furthermore let 

KL = KL (LW) = Ke -^ + 0(^) . (2.44) 

In addition, parametrising the leading L-dependence of the screening masses via 

/4 = J2 ' alld = J2 ' ( 2 - 45 ) 

we obtain from the leading corrections in (2.41) and (2.42), respectively, 

Iz( K c) /„ -i \ K c(2 + « c ) /ri AC\ 

Ck 77"^ = 2/tc ~ 1 Cz ~ 771 7 Cg ' 2 - 46 

7z(«c) 6(1 - k c ) 

/^(« c ) _ (2« c -1)(2 + k c ) / k\ (2 + Kc ) \ 

° K / G («c) 4(1 - Kc ) CZ ^ c (1 + Kc )(l - « C )J CG • 

The fact that we have 3 unknowns, c K , Cz and Cq, from only 2 equations here reflects an 
ambiguity in our parametrisations (2.20) of the finite-volume infrared behaviour of the gluon 
and ghost dressing functions in the momentum regime of (1.1). As a parametrisation of the 
leading corrections to the infinite- volume infrared forms with originally 2 parameters in each of 
the dressing functions, Zi(k 2 ) and Gl(A; 2 ), Eqs. (2.20) are somewhat redundant, of course. Even 
after eliminating one of the originally 2 exponents, the residual redundancy is manifest here. 
An L 2 k 2 independent exponent is excluded, however, as c K = implies that also cz = c G = 
for k c = (93 - VT201) /98. 

In fact, in units of c K , Eqs. (2.46) yield numerical values 

c z ~ 238 c K , c G ^86c K . (2.47) 



Because c K > for positive c z / G (c./., Eqs. (2.45)), with the definition of c K in (2.44), we find 
that Kl approaches k c from below. 

Moreover, Eqs. (2.47) imply that the momentum dependence of the exponent KL(L 2 k 2 ) is a 
rather weak effect in comparison to the influence of the screening masses, which thus account 
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for the dominant finite- volume corrections at sufficiently large k 2 L 2 . Therefore, a quantitatively 
good approximation for large volumes and momenta in (1.1) is given by 



Z L (k 2 ) = Z m 



Aqcd 
k 2 + » 2 z , 



k 2 



k 2 + IM% 



2 ' 



G L {k 2 ) = Gj A ^ CD 



k 2 + ^ G , 



(2.48) 



with the infinite volume k's. With this approximation of neglecting the momentum dependence 
of the infrared exponent kl, estimates of the screening masses Hz/g can m principle be obtained 
from Eqs. (2.31) and (2.34), (or (2.39)) at k 2 = rather than considering large k 2 for the 
leading l/(L 2 k 2 ) corrections as done so far. Because the resulting equations will then necessarily 
involve the summations over the discrete loop momenta, such an estimate would already need 
a numerical solution to these equations. Here, we rather present our full numerical solutions to 
the torus DSEs in the next section. The masses extracted from these full solutions in Sec. 3.4 
are in good agreement with the ratio cz/cg determined here. 

Finally, inserting the screening masses Hz/g — c z/g/L 2 with L-independent coefficients cz/g 
into the IR-asymptotics (2.48), for vanishing momentum k 2 with 5Z = 5G = and —kz l /2 = 
k Gl ps k c in (2.20), we obtain an estimated behaviour with L of the infrared limit of the gluon 
and ghost dressing functions as 



lim 



Z L (k' 



fc 2 ^o k 2 



oc L 2 * 1 - 2 *") , 



lim G L (0) oc L 

fc 2 ^0 



2h r 



[2A9) 
[2.50) 



The L-dependence of the dressing functions at k — mimics the infrared leading momentum 
behaviour. In particular, as (1 — 2k c ) ps —0.19, this implies that the zero momentum gluon 
propagator is expected to decrease extremely slowly with the volume V = L A towards zero in 
the infinite volume limit oc V~ e with an exponent e = (2k c — l)/2 which is smaller than 0.1. 
The infrared limit of the ghost propagator, on the other hand, increases faster than linear with 
the length L. We will compare these qualitative features to the full numerical solution in the 
next section. 
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Fig. 2. Sketch of the momentum grid dual to the four-torus. The summation over complete hyper- 
spheres is indicated by fully drawn circles. The hyperspheres depicted by dashed lines are not complete 
if one uses Cartesian cutoffs instead of an 0(4) invariant one. 

3 Numerical results 

3.1 The momentum lattice 

As already mentioned in the previous section, in a 4-dimensional hypercubic volume V = L 4 
with periodic boundary conditions, the momentum integrals of the infinite volume DSEs are 
replaced by sums, c.f., (2.6). For the numerical treatment of the equations it is convenient to re- 
arrange these summations in a spherical coordinate system [17]. A two-dimensional illustration 
of this procedure is given in Fig. 2. We thus write (2.6) in the form 

^ E (■■■) = (■■■), ( 3 - 51 ) 

n e Z 4 rn, 3 

where m labels sums over hyperspheres, each containing all momentum vectors of the same 
absolute value and j numbers the individual vectors on a given hypersphere. The resulting 
double sum corresponds to the splitting between radial and angular integrals in the infinite 
volume DSEs. This correspondence is very good for large hyperspheres (i.e., in the ultraviolet 
momentum region), where a large number of vectors on a given sphere samples the angular 
dependence of the integral kernels well. However, in the infrared, i.e., on the innermost spheres, 
this sampling will be poor, thus resulting in hypercubic artefacts. We will identify these in our 
solutions in Sec. 3.3. Note, however, that the sole dependence of the dressing functions on the 
squared momenta is not touched by these artefacts. We explicitly verified that the dressing 
functions Z{p 2 ) and G{p 2 ) do not depend on the direction of the momentum p. 

For numerical reasons we also have to introduce a momentum cut-off that limits the extent of the 
momentum lattice. In the infinite volume formulation such a cut-off Auv is 0(4)-invariant. The 
resulting renormalised dressing functions are independent of the value of Auv, as they should be 
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in a renormalisable quantum field theory. In general, this is no longer true on a four-torus and we 
expect cut-off artefacts. It turns out that these effects are sizeable only when 0(4)-symmetry is 
badly broken as is the case, e.g., when restricting the sums in (2.6) in each Cartesian direction. 
Much better results are obtained, if an (approximately) 0(4) -invariant cut-off is introduced by 
using the summation procedure (3.51) and summing only over such hyperspheres which do not 
receive any further vectors when enlarging the momentum grid. In Fig. 2 these are indicated by 
the solid lines as opposed to the dashed ones. Since these fully occupied outer spheres typically 
contain a large number of momentum vectors, the corresponding 'angular sum' over j is close 
in value to the corresponding angular integrals in the infinite volume limit and is thus a good 
approximate representation of the corresponding 0(4)-symmetry. Note that we have used a 
similar procedure in our infrared analysis around and below Eq. (2.38). 

3.2 Renormalisation procedure 

Further details of our numerical method to solve the coupled system of Eqs. (2.8), (2.7) have 
been described in [17,18,35] and shall not be repeated here. The only major difference to 
previous calculations of DSEs on the torus is a modification of the renormalisation conditions. 
These modifications turn out to have a substantial impact on the large volume limit. 

The propagator DSEs (2.10) and (2.11) resulting from the renormalisation condition (2.9) 
in Sec. 2 make the approach towards the infinite volume solutions at large torus volumes 
transparent. A direct numerical implementation of eqs. (2.10) and (2.11) requires a sufficiently 
large UV-cutoff A uv . Unfortunately, for reasons of CPU-time we are restricted to a rather low 
UV-cutoff of the order of Auv = 2 — 3 GeV. For such a low cutoff one encounters convergence 
problems in the infinite volume DSEs which prevent a reliable extraction of the ghost and gluon 
renormalisation factors Z 3 and Z 3 . We therefore adjust the renormalisation procedure of Sec. 2 
to the numerics as explained in the following. 

We substitute the infinite volume subtractions in eqs. (2.10) and (2.11) by subtractions on the 
torus at a fixed scale s 2 , which in general does not have to be equal to the renormalisation 
point. In the symbolical notation 



1 



Z 3 + g 2 N c U ghost (p 2 ) , 



(3.52) 



G( P 2 ) 



1 



Z 3 + g 2 N c U glue (p 2 ) , 



(3.53) 



Z(p 2 ) 



for the equations (2.7) and (2.8) this procedure yields 




(3.54) 
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l^r + g 2 N c (u glue (p 2 ) ~ IW(s 2 )) , (3.55) 



Z(p 2 ) Z(s 2 ) 

similarly to (2.10) and (2.11). The renormalisation constants Z 3 and Z 3 are traded for two 
input values at G(s 2 ) and Z(s 2 ). These values are not independent of each other. In the infinite 
volume limit this can be seen from the Slavnov- Taylor identity 

Z^Alv) = Z g ^\kl v )Z^\kl w )Zl /2 ^\kl Y ), (3.56) 



In Landau gauge, the ghost-gluon vertex is ultraviolet finite, and therefore one always has the 
choice of setting Z 1 — 1 [26]. Thus for a fixed Z g , corresponding to a fixed choice of the renor- 
malised coupling a(/i 2 ) = g 2 (fi 2 ) / An , one has a unique relation between the renormalisation 
factors Zj, and Z 3 of the ghost and gluon propagators. In terms of renormalisation conditions 
this corresponds to the relation 

G^fZ^ 2 ) = 1, (3.57) 



i.e. the ghost and gluon dressing functions may not be renormalised independently for a given 
renormalised coupling [25]. Accordingly, the input values G(s 2 ) and Z(s 2 ) have to be chosen 
to satisfy (3.57). In the infrared analysis in Sec. 2 these could be taken from the well-known 
infinite volume infrared solutions in accordance with the renormalisation conditions (2.9). For 
a numerical solution, however, this is not an option. 

The simplest way to implement (3.57) would be to choose s 2 = /i 2 , fix Z(s 2 ) and G(s 2 ) = 
l/Z(s 2 ) 1 / 2 . Unfortunately this is not possible: the function g — > /i 2 is not known before the 
equations are solved completely and accordingly choosing a specific g is not sufficient to de- 
termine the value of /x 2 . One thus chooses an arbitrary, usually large, subtraction point s 2 , 
fixes Z(s 2 ) and then solves for a range of values for G(s 2 ). In the infinite volume limit one 
obtains continuous and differentiable solutions only for an extremely narrow range of values 
G(s 2 ) around the correct one implementing (3.57) [34]. 

However, this is no longer true on the compact manifold. We found that for given values of g(fi 2 ) 
and Z(s 2 ) one can generate a continuous array of solutions by varying G(s 2 ). These solutions 
all behave different in the infinite volume limit. In previous works [17,18] these ambiguities have 
been resolved by reading off Z(s 2 ) and G(s 2 ) at a large s 2 from the infinite volume solution. 
However, this procedure is not sufficient. Since one works with a fixed cutoff Auv o n the torus, 
one encounters 0(1/Auv) effects in the ultraviolet momentum region, similar to 0(a)-effects in 
lattice QCD. Thus the ultraviolet behaviour of the torus solutions is slightly different from the 
infinite volume results. This difference is significant for the renormalisation procedure and in 
turn affects the infrared behaviour of the solutions and also the scaling behaviour with volume. 

To avoid these UV-cutoff effects we modify the RG procedure by utilising the analytic results in 
Sec. 2: the RG condition on the torus is adjusted to approach the infinite volume RG condition 
in the limit of very large volumes as implemented in Sec. 2 in (2.10) and (2.11). To that end we 
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choose G(s 2 ) for fixed g and Z(s 2 ) accordingly. As a result we find that the values of Z(s 2 ) and 
G(s 2 ) are different but still close to the corresponding values of the infinite volume solution. 

To summarise: we find that 0(a)-effects have to be accounted for in the renormalisation proce- 
dure for DSEs on the torus to obtain the correct infinite volume limit. This technical improve- 
ment compared to previous works [17,18] allows us in turn to reliably explore the large volume 
behaviour of the ghost and gluon propagators. 

3.3 Propagators 

Our numerical results for the ghost and gluon propagator on different volumes are shown in 
the two graphs of Figure 3. The momentum scale is fixed by comparison with corresponding 
lattice calculations: we demand that the maximum of the gluon dressing function in the DSE- 
approach occurs at the same momentum scale as in the lattice dressing function. Hereby we 
assume implicitly that truncation errors due to the neglect of the gluonic two-loop diagrams 
only mildly affect the position of the maximum. As a crosscheck we compare the resulting 
values of the DSE and the lattice gluon propagator at the lowest momentum point accessible 
on manifolds of similar volume. The result, shown explicitly in Sec. 3.5 below, suggests that 
scale uncertainties relative to the lattice scale might roughly be of the order of 10 percent. 

We discuss results on seven different volumes V = L 4 ; the corresponding box lengths L are given 
in the legends of Figure 3. One clearly observes that the infinite volume solutions of the gluon 
propagator D(p 2 ) and ghost dressing function G(p 2 ) are more and more approached by the torus 
solutions with increasing volume. This tendency will be quantified in Sec. 3.4. Qualitatively one 
can see the following behaviour from the plot: the propagator seems to be divergent at volumes 
of V ?s (4 — 8 fm) 4 . Note that these volumes are already large compared to the ones used in most 
lattice calculations of observables. For even larger volumes the propagator bends downwards to 
reach a plateau at roughly V ~ (9 fm) 4 . For volumes larger than V ~ (10 fm) 4 the propagator is 
infrared vanishing and therefore qualitatively similar to the infinite volume limit given in (2.2) 
and (2.3). 

This finding becomes even more pronounced when one applies momentum cuts: from the inset 
of the graph for the gluon propagator one can see that the first three momentum points on 
the dual torus behave differently than the other points. In the language of Sec. 3.1 this means 
that the 'angular sum' over j on these spheres is by no means a good representation of the 
corresponding angular integrals in the infinite volume limit. If one omits these points in the 
representation of D(p 2 ) = Z(p 2 )/p 2 and G(p 2 ), then the turnover of the gluon propagator shifts 
from roughly V ~ (9fm) 4 to V ~ (11 fm) 4 . In such a representation only the propagator of 
our largest volume, V = (13.8 fm) 4 , can be seen as infrared vanishing. Taken at face value this 
means that contemporary lattice calculations are still far away from the critical volume where 
the behaviour of the infinite volume solution can be deduced. 

Similar conclusions can be drawn from the results for the ghost propagator. In Figure 3 we 
observe that the first two or three points on each result for the ghost bend away from the 
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Fig. 3. Numerical solutions on tori with different volumes compared to the infinite volume limit. The 
upper graph shows the gluon propagator, whereas on the lower graph the ghost dressing function is 
depicted. 
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power law behaviour of the infinite volume solution. The remaining curves have more and more 
points on the 'scaling region' of the ghost dressing function, where the infinite volume power 
law develops. The critical volume, where the full infinite volume power law can be seen on the 
torus, is of roughly the same size as the one for the gluon propagator. 

It is interesting to note that the lowest momentum point of the gluon and ghost solutions follow 
a curve distinct from the infinite volume solution. It turns out that in both cases this curve can 
be described by the infinite volume solution multiplied by a constant factor (the dotted lines 
in the graphs). A similar observation can be made for the second lowest point. As yet we have 
no full analytical understanding of this observation. 

This behaviour has interesting consequences for the running coupling, defined in eq. (2.4) and 
shown in Figure 4: at least the first two points of every solution on a torus will always remain 
significantly below the infinite volume solution no matter how large the volume is. Thus on every 
compact manifold one finds a running coupling which looks infrared suppressed, although the 
corresponding coupling from the infinite volume limit has an infrared fixed point. Taken at 
face value this result means that it is extremely difficult if not impossible to verify infrared 
fixed points of the running coupling in lattice calculations. Indeed, all recent determinations of 
vertex couplings from lattice QCD find infrared vanishing [36] or strongly suppressed couplings 
[37] whereas the infinite volume analysis unambiguously predicts fixed points [24,1,31,32]. Our 
observation may well serve to explain this discrepancy. 

Finally, in the lower panel of Figure 5 we present results for two tori with the same physical 
volume but different ultraviolet cut-offs given in the legend. The two corresponding solutions for 
the gluon propagator are close to each other, although one observes a slight distortion of shape 
from the curve with the small cut-off compared to the one with the larger cut-off and the infinite 
volume solution. From this result one may conclude that cut-off effects in the DSE-solutions on 
a torus are small and do not significantly affect the infrared behaviour of the propagators^ 

3-4 Emergence of the infrared asymptotic behaviour 

We now study the rate of approach towards the infinite volume solution in some more detail. 
Because of the redundancy in the parametrisations (2.20) as discussed in Sec. 2.4, direct fits 
to our data simultaneously of all free parameters in the finite-volume infrared forms are not 
possible. Fits of both screening-mass parameters would in principle be possible from the expres- 
sions in Eqs. (2.48). It will turn out, however, that the volumes in our numerical results which 
range up to approximately (14 fm) 4 are nevertheless still too small for the leading infrared 
exponent k c , or any other momentum and volume independent constant, to describe the data 
in a reasonable range of momenta. For these volumes the available infrared momenta all lead 
to values of l/(k 2 L 2 ) which are beyond the range of validity of the expressions in (2.48), where 
the exponent k c is already close to its infinite volume value but the screening masses fiz/G are 

4 In the light of these results we also conclude that the larger cut-off effects noted in [18] are artefacts 
of the renormalisation procedure used therein. 
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Fig. 4. The running coupling a(p 2 ) = a(fi 2 )G 2 (p 2 )Z(p 2 ) as defined from the ghost-gluon vertex. 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



— Infinite volume 

— A=2.2 GeV (L = 13.5 fm) 

— A=1.7 GeV (L = 13.2 fm) 




Fig. 5. The gluon propagator calculated on tori with the same physical volume but different ultraviolet 
cut-offs. 

still present. We will find, in fact, that there is not yet a clear separation of scales (1.1) for any 
momentum in these volumes in the first place. 
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We therefore adopt two different models, each to fit the data in two steps. To illustrate this 
procedure we start from the ghost dressing function for which it yields the more stable results 
in the available volumes than for that of the gluon. Instead of (2.48) we assume a ghost dressing 
function at low momenta of the form, 

(A2 \ Aghast 

wfk) • (3 ' 58) 

in which we allow both, the screening mass /iq and the effective exponent K g h os t to depend on 
the volume, i.e., on L. We then proceed in the following two alternative ways: 

(1) In this model we first set the screening mass to zero, i.e., \iq = in (3.58), and determine 
the effective exponent K g h ost by fitting the third to fifth lowest momentum values of the 
solution for every given volume to the resulting pure power law oc (AQ CD /A; 2 ) K9 ' iOSt . The 
results are plotted over the inverse box length as the triangles in the left panel of Fig. 6. 

We then keep Kghost fixed and fit in a second step the remaining parameters on a larger 
region. Here we typically use the lowest eight momenta. The resulting values for ghost 
screening mass are the triangles in the left panel of Fig. 7 which all lie on a straight line 
cx 1/L with offset zero. 

The resulting fits are stable with respect to moderate variations of the fit regions in each step 
(within error bars) and nicely reproduce the same results. Note that this method predominantly 
absorbs the finite volume effects in the effective exponent. The screening mass essentially ac- 
counts only for the remaining mismatch after this first step. 

(2) In this fit-model we adopt the extremely opposite point of view. We first obtain an effective 
exponent K g h ost which is completely independent of any finite-volume effects by fitting the 
pure power law form to the infinite volume solution in the same momentum range used in 
the first step of model (1) above, that of the third to the fifth lowest momentum values of 
a given torus solution. The resulting values are plotted for comparison as the dots which 
lie above the triangles, also in the left panel of Fig. 6. 

Fitting the remaining parameters to that torus solution then essentially works as in model 
(1) above. We use the same larger momentum region of the typically eight lowest values, 
and the fits are similarly stable w.r.t. variations of that region as above. The resulting 
screening masses are larger than in model (1) and shown as the squares in Fig. 7. 

The crucial distinction is that in this second model the finite-volume effects are entirely ac- 
counted for by the screening mass alone. For the volumes available the effective exponent n ghost 
in (3.58) is not equal to k c yet. It changes with the volume because the momentum range where 
it is obtained does. Nevertheless, this effective exponent is obtained from the infinite volume 
solution and thus independent of finite volume effects. It rather describes the deviations from 
the leading infrared behaviour already of the infinite volume solution in the corresponding mo- 
mentum range, due to the fact that these momenta are not yet far enough below Aqcd- The 
volumes are too small to fit a range of momenta (1.1) in. 
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Fig. 6. The infrared exponents K g h os t an d Kgiue from fits to the torus solutions on different volumes, 
c./., model (1) as described in the text, compared to the effective power law fits to the infinite volume 
solution in the same momentum region for model (2). 

The fact that both alternative fit-models are possible reflects the redundancy in our finite- 
volume infrared parametrisations as mentioned above. In volumes which are not quite large 
enough for the leading volume asymptotics we can neither describe the data with constant 
exponent nor without screening mass. Both, an effective exponent and a screening mass are 
necessary in these volumes. However, to allow this both underconstrains the fits. The two 
extreme models to fit the data by absorbing a maximum of finite-volume effects either in the 
effective exponent or the screening mass serve as a measure of this redundancy. 



The results approach each other and this redundancy decreases with increasing volumes as 
expected. In the volumes of the order of up to approximately (14 fm) 4 presented here the 
mismatch between both models does not quite vanish yet. We can extrapolate the data from 
both fit-models to larger volumes, however, to make a rough prediction about what volumes 
are necessary for the unique leading finite-volume infrared behaviour as parametrised in (2.48) 
to set in. The result gives the following quite consistent overall picture: 

Comparing Figs. 6 and 7 we observe that for box lengths of approximately L = 40 fm (i.e., 
1/L « 0.025 fm -1 ) the effective exponents from torus and infinite volume solutions approach 
each other and eventually meet extremely close to the asymptotic infinite-volume value k c « 
0.595. For lengths above 40 fm we therefore expect Kghost — k c to be a very good approximation. 
The screening masses are not zero yet, but also approach each other very well at volumes of this 
size. This means that from our extrapolations we would estimate that the asymptotic forms 
(2.48) can be expected to describe the finite volume infrared behaviour well for volumes of L = 
40 fm in size and more. 

As described in the previous section, the dressing function of the gluon propagator is more 
difficult to describe by such fits in the intermediate volumes here because its momentum de- 
pendence deviates much more strongly from its infinite-volume infrared asymptotics even at 
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the lowest momenta available in our volumes. For the gluon dressing function model (1) only 
leads to stable results for the screening mass \iz in our largest volume V = (13.8 fm) 4 which is 
included in the right panel of Fig. 7. The volume dependence of the gluon dressing function is 
stronger than that of the ghost. This is apparent from the right panel in Fig. 6 for model (1), 
with the volume dependence predominantly accounted for by the effective exponent which is 
obtained from the fiz = fits of the form, 

Z L (k 2 ) = Z m -T2— , (3.59) 

\ A QCD/ 



and which depends more strongly on 1/L than it does for the ghost. But it again approaches 
the effective exponent of the infinite volume solution at box lengths of about 40 fm. 

In model (2), with no finite- volume effects in the effective exponents, we observe the same 
stronger volume dependence in the gluonic screening mass /iz as compared to hq m Fig- 7. 
Again, however, within the considerable errorsjjj] the extrapolations of the data for the gluonic 
screening masses from models (1) and (2) are consistent with an approach of the two at about 
the same length scale of L = 40 fm. 

Note that the screening- masses of model (1) are extrapolated using linear fits (of the fits) 
inversely proportional to the length, 

li z/G cx 1/L . (3.60) 



In model (2) on the other hand, our extrapolations of the screening masses are exponential in 
nature from fits of the form 

Hz/a oc m ( exp{// L} - 1) , (3.61) 

with dimensionful constants m (in MeV) and I (in fm). In our intermediate volumes both forms 
differ considerably. At sufficiently large L (with our estimate > 40 fm), however, they both 
agree with the scaling Ansatz (2.45) used the infrared analysis of Sec. 2.4. 

As a final cross-check, for the ratio of the screening masses, 

/i z //i G « v/238/86 « 1.7 , (3.62) 



as obtained from Eqs. (2.47) in Sec. 2.4 for asymptotically large volumes, here we roughly 
obtain fiz/^G ~ 1-6 for their ratio from the corresponding ratio of the slope of the linear fits in 
model (1), the solid lines through the triangles in Fig. 7. This is in perfect agreement with the 

5 For the gluonic screening-mass in model (1) we simply draw a straight line between our only value 
at L = 13.8 fm and zero at 1/L = in Fig. 7 (right). 
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P min [MeV] P min [MeV] 
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Fig. 7. In the left panel the behaviour of the ghost mass [iq with the inverse box length 1/L is shown. 
The error of five percent reflects the systematic uncertainty of our fit procedure. The straight line is 
a linear fit with zero offset. In the right panel the same is shown for the gluon mass fiz- However, the 
systematic uncertainty is larger, about forty percent, see text. 

predicted asymptotic ratio (3.62), but the procedure has a large error of around 50%, mainly 
from the single value [iz = 70(30) MeV for the (13.8 fm) 4 volume which determines the slope 
in the gluonic screening mass in model (1). 

For the effective exponents of both dressing functions we always use extrapolations (the long 
dashed lines for model (1) and the short dashed lines for model (2) from the infinite volume 
solutions) based on quadratic polynomials^! in 1/L with the 1/L = value fixed at k c « 0.595, 

K fit (l/L) = k c - -j- - . (3.63) 



For the dressing functions in model (1) we have also included the constant term Kfu{0) in the 
quadratic fit of the form (3.63): From the ghost dressing function this yields K/«t(0) = 0.597(3) 
in perfect agreement with the analytic result of k c = 0.595 at 3-digit precision. 

Note that for the gluon dressing function, the model (1) fits using the form (3.63) with fixed 
k c produce a coefficient a\ of the linear term consistent with zero. Including the constant term 
in the complete 3 parameter fit is then not stable anymore. It requires dropping the linear 
term, i.e. freezing a\ = 0. Doing so, however, one then obtains Kfu(0) = 0.60(2) which is also 
consistent with the analytic result. 

Finite volume corrections to k c of the form (2.44) used for our infrared analysis in Sec. 2.4 can 
not reasonably be extracted from our present data. This is not unexpected because it requires 

6 The only exception here is the effective exponent of the infinite-volume gluon solution, the short 
dashed line in the right panel of Fig. 6, which requires the inclusion of a cubic term. 
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volumes of a size at which the scaling relation (2.29), i.e., k,z l + ^ k g l — 0, for the infrared 
exponents is already valid. Here, c.f., the definition of the effective exponent for our fits in 
(3.59), this is when 

l^glue ~ foghost j 

(3.64) 

which according to our estimate again requires L pa 40 fm, c.f., Fig. 6, and is thus still relatively 
far beyond our numerical data to attempt reasonable fits and extrapolations. 

We furthermore note that the length scale where the behaviour of the gluon propagator turns 
from k < 0.5 (infrared diverging) to k > 0.5 (infrared vanishing) is of the order of at least 
L = 6 fm for model (2) (from the infinite volume effective exponent) up to L = 10 fm for model 
(1). The exact turning point might also be somewhat truncation dependent. 

In summary, we have seen in the previous two sections that the qualitative infrared behaviour 
of the gluon and ghost propagators of Landau gauge QCD starts to emerge in volumes of 
estimate sizes between L = 10 fm and 15 fm. To reliably extract infrared exponents and other 
quantitative results about their infrared behaviour, one furthermore needs a certain range of 
momenta with a clear separation of scales (1.1) which requires much larger volumes. The leading 
finite-size effects are under good control, and can be described by the simple screening masses 
\xzig oc 1/1/ via Eqs. (2.48), when the effective exponents approach their infinite volume scaling 
behaviour and are sufficiently close to k c . Our extrapolation method predicts that this might 
require volumes of about 40 fm in size. 

3.5 Comparison to lattice results 

We now compare our results from DSEs on a torus to the ones from lattice calculations. Lattice 
studies of the Landau gauge gluon propagator have quite a long history already [38,39,40,41,42]. 
Nowadays, lattice data for the gluon propagator is available on impressively large lattices. 
The authors of [10] report on an S77(2)-study on a 52 4 -lattice, whereas in [11] results from 
an SU(3) calculation on a 56 4 -lattice are discussed with larger lattices on their wayf^l The 
physical volumes of both these studies are comparable. The resulting propagators are very 
similar with small deviations for the first few points in the infraredjjj In the upper graph of 
Figure 8 we display the SU (3)-results together with data on a smaller volume [45] and compare 
with DSE-results on tori with similar volumes. 

The qualitative agreement of the solutions at similar volumes in the infrared is striking]^] 

7 For results from simulations using large asymmetric lattices, see [43]. 

8 This independence of the number of colours has also been found in an earlier study [44]. However, 
there a different gauge, the LaPlacian gauge, has been used. 

9 Differences in the intermediate momentum regime around 1 GeV are truncation artefacts of the 
DSE solutions. One can show analytically that this is exactly the region where the omitted gluonic 
two-loop diagrams contribute significantly. 
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Fig. 8. DSE results for the gluon propagator and ghost dressing function on tori with different vol- 
umes compared to recent lattice calculations on similar manifolds. The lattice data are taken from 
Refs. [11,45]. 

Whereas both, the lattice and the DSE result at the smaller volume V ~ (4.6 fm) 4 seem to 
diverge, one starts to observe an infrared finite one, or perhaps even a slight infrared suppression, 
at the larger volumes V ~ (9.7 fm) 4 . This indicates that the scaling behaviour of the lattice 
results with volume may be very similar to the ones of the DSE solution. If so, then our results 
from the last section predict also a turnover of the lattice data when even larger volumes are 
considered. 
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The situation is less clear for the ghost dressing function, which was first studied on the lattice 
in [46]. Our results for three different volumes, V = (3.9, 4.6, 9.7 fm) 4 , are compared to the 



SU(3) lattice results of [11.45]] 10 1 For the DSE solutions we observe a characteristic deviation 
of the two lowest momentum points at each volume from the infinite volume solution. As shown 
in the previous section these deviations correspond to a ghost mass which goes to zero in the 
infinite volume limit. The lattice results do not seem to show such behaviour as yet. Even 
though the lattice volumes herein are roughly between 3 fm and 4.5 fm, and thus still rather 
small compared to our analysis in the previous section, there appears to be not much sign of a 
volume dependence at all at this point. The situation is indeed reminiscent of the corresponding 
torus solution from the previous DSE studies. Perhaps revisiting the renormalisation procedure 
as described in Sec. 3.2 might have a similar influence in the analysis of the lattice data. 

In these relatively small volumes, however, some finite- volume effects in our DSE solutions might 
also still be somewhat truncation dependent and could therefore contribute to the differences 
between the DSE solutions and on the lattice simulations. Furthermore, effects from Gribov 
copies are known to influence the lattice results for the ghost propagator much stronger than for 
the glue. This has been investigated by comparing results from gauge fixing to the first Gribov 
copy (fc) to the results of a more involved procedure of repeated gauge fixing and selecting 
then the best copy (be). On rather small volumes these two procedures have been compared 
in [45] and noticeable effects observed. Similar conclusions have been drawn in [48]. Therefore, 
a meaningful comparison between the DSE and lattice ghost propagator not only awaits an 
improved truncation scheme in the DSEs and results on larger volumes on the lattice, but also 
a further clarification of gauge fixing effects in both, the continuum and lattice studies. 



4 Conclusions 



The results presented here explicitly demonstrate that Landau gauge propagators from lattice 
and continuum calculations are compatible, provided finite volume effects are properly taken 
into account. The approach of the torus solutions towards the infinite volume limit is slow, 
but steady and quantifiable. Our results show that the characteristic maximum in the gluon 
propagator should be measurable in volumes of the order V ~ (10 — 15 fm) 4 which come into 
grasp for lattice calculations. This important qualitative effect will be the hallmark of the onset 
of the asymptotic infrared region. 

The results also demonstrate, however, that extremely large volumes will be necessary for a 
more quantitative comparison between results in the continuum and on the lattice. The reason 
is, of course, the requirement (1.1) to observe the asymptotic infrared behaviour. Hence, even 
verifying the scaling of the infrared exponent according to figure 6 requires spanning an enor- 
mous range of volumes. Thus the agreement between both methods can be at best qualitative 
for quite some time to come. That so far not even a very good qualitative agreement could be 

10 For corresponding results in 577(2) see [47]. To date no systematic lattice study on the effects of 
the gauge group on the ghost propagator exists. 
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observed can also be understood with the results presented here. The most characteristic quali- 
tative feature of the gluon propagator is its maximum. To observe it, it is trivially necessary to 
also observe its infrared suppression. This suppression depends on the quantitative value of the 
infrared exponent K g i ue - In particular, its value in Landau gauge, compared, e. g., to Coulomb 
gauge, leads to a very weak infrared suppression. Hence very large volumes are necessary to 
observe this suppression, and hence the maximum. 

Nonetheless, if a clear maximum emerges in the gluon propagator on larger lattices, confirming 
the predictions made here, it will definitely be an important piece in our understanding of the 
confinement mechanism in Landau gauge. 

Finally we would like to remark that all volume effects discussed here are only typical for 
correlation functions involving intrinsically massless fields. Objects, which have an intrinsic 
mass-scale, like quarks or hadrons, will only indirectly be affected by such finite volume effects. 
For quarks these effects may be still significant, as can be seen from the results of [19]. The 
situation is different, however, for hadrons: the specific far infrared behaviour of the gluon 
correlation functions are expected to be nearly irrelevant to hadronic observables [1,23,49]. 
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